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We define a new measure of quantum correlations in bipartite quantum systems given by tlie 
Bures distance of tlie system state to tlie set of classical states witli respect to one subsystem, 
that is, to the states with zero quantum discord. Our measure is a geometrical version of the 
^-^ . quantum discord. As the latter it quantifies the degree of non-classicality in the system. For pure 

^vj ' states it is identical to the geometric measure of entanglement. We show that for mixed states it 

coincides with the optimal success probability of an ambiguous quantum state discrimination task. 
Moreover, the closest zero-discord states to a state p are obtained in terms of the corresponding 
optimal measurements. By solving the quantum state discrimination problem for two qubits we 
determine the geometric quantum discord and the closest zero-discord states explicitly for states p 
with maximally disordered marginals. 



m 



Oh: 

< 



Ph. 



> 

m 
m 
m 

o 
en 



X 



PACS numbers: 03.67.Mn, 03.67.Hk 



Dominique.Spehner@ujf-grenoble.fr 



I. INTRODUCTION 

One of the basic question in quantum information theory is to understand how quantum correlations in composite 
quantum systems can be used to perform tasks that cannot be performed classically, or that lead classically to 
much lower classical efficiencies [1]. These correlations have been long thought to come solely from the entanglement 
among the different subsystems. This is the case for quantum computation and communication protocols using pure 
states. For composite systems in mixed states, however, there is now increasing evidence that other types of quantum 
correlations, captured by the quantum discord of Ollivier and Zurek [2] and Henderson and Vedral [3], could provide 
the main resource to exploit in order to outperform classical algorithms [4-6] or in some communication protocols [7- 
9]. The quantum discord quantifies the amount of mutual information not accessible by local measurement on one 
subsystem. One can generate mixed states with nonzero discord but no entanglement by preparing locally statistical 
mixtures of nonorthogonal states, which cannot be perfectly discriminated by measurements. The strongest hint 
suggesting that the quantum discord may in certain cases quantify the resource responsible of quantum speedups is 
provided by the deterministic quantum computation with one qubit (DQCl) of Knill and Laflamme [10] leading to 
exponential speedups with respect to known classical algorithms. Even though no entanglement is present during the 
whole computation ]11], a nonzero quantum discord appears between the pure single qubit and the other qubits [4, 5]. 
This can be nicely interpreted by using the monogamy relation ]12] between the quantum discord of a bipartite system 
AB and the entanglement of B with its environment E if ABE is in a pure state ]6]. 

A mathematically appealing way to quantify quantum correlations in multipartite systems is given by the minimal 
distance of the system state to a separable state [13]. The Bures metric ]14, 15] provides a nice distance cIb on the 
convex cone of density matrices, which has better properties than the Hilbert-Schmidt distance ^2 from a quantum 
information perspective. In particular, ds is monotonous and riemannian [16] and its metric coincides with the 
quantum Fisher information [17] playing an important role in high precision interferometry. As a consequence, the 
minimal Bures distance to separable states satisfies all criteria of an entanglement measure [13], which is not the 
case for the distance d2. This entanglement measure has been widely studied in the literature ]18-20]. By analogy 
with entanglement, a geometric measure of quantum discord has been defined by Dakic, Vedral, and Brukner [21] as 
the minimal distance of the system state to the set of zero-discord states. This geometric quantum discord has been 
evaluated explicitly for two qubits ]21]. However, the aforementioned authors use the Hilbert-Schmidt distance ^2- 

The aim of this work is to study a similar geometric quantum discord (GQD) as in Ref. ]21] but based on the 
Bures distance rfs, which seems a more natural choice. This distance measure of quantum correlations has moreover 
a clearer geometrical interpretation than other measures ]13, 22, 23] based on the relative entropy, which is not a 
distance on the set of density matrices. We show that it shares many of the properties of the quantum discord. Most 
importantly, our geometrical approach provides further information not contained in the quantum discord itself. In 
fact, one can look for the closest state(s) with zero discord to a given state p, and hence learn something about the 
"position" of p with respect to the set of zero-discord states. We calculate the Bures-GQD explicitly for pure states 
of bipartite quantum systems and for mixed two-qubit states with maximally mixed marginals. In both situations 
we also determine the closest zero-discord state(s) and compare our results to those obtained by using the Hilbert- 
Schmidt distance. Our results are obtained by first establishing a link, which holds for arbitrary mixed states p in 
any bipartite system, between the Bures-GQD and the closest zero-discord states and the well-known minimal-error 
quantum state discrimination problem. 

The task of discriminating states pertaining to a known set {pi, • • • , p„} of states pi with prior probabilities rji 
plays an important role in quantum communication and quantum cryptography. For instance, the set {pi, • • • ,/0„} 
can encode a message to be sent to a receiver. The sender chooses at random some states among the p^'s and gives 
them to the receiver, who is required to identify them and henceforth to decode the message. With this goal, the 
receiver performs a measurement on each state given to him by the sender. If the pi are non-orthogonal, a perfect 
discrimination is impossible, so that the amount of sent information is smaller than in the case of orthogonal states. 
The best the receiver can do is to find the measurement that minimizes in some way his probability of equivocation. 
Two distinct strategies have been widely studied in the literature (see e.g. the review ]24]). In the first one, the 
receiver seeks for a generalized measurement with {n + I) outcomes, allowing him to identify perfectly each state pi 
but such that one of the outcomes leads to an inconclusive result (unambiguous quantum state discrimination). The 
probability of the inconclusive outcome must be minimized. In the second one, the receiver looks for a measurement 
with n outcomes yielding the maximal success probability Ps = "^i^iVi-Piliy where P^ij is the probability of the 
measurement outcome i given that the state is pi. This strategy is called minimal error (or ambiguous) quantum 
state discrimination. The maximal success probability and optimal measurements are known explicitly for n — 2 ]25], 
but no general solution has been found so far for more than two states (see, however, [26]) except when the pi are 
related to each other by some symmetry and have equal probabilities rji (see [24, 27] and references therein). However, 
several upper bounds on Pg^ are known ]28] and the discrimination task can be solved efficiently numerically ]29]. 
Let us also stress that unambiguous and ambiguous discrimination have been implemented experimentally for pure 



states [30] and, more recently, for mixed states [31], by using polarized light. 

Our main result shows that the Bures-GQD of a mixed state p is equal to the maximal success probability Pg in 
ambiguous state discrimination of a family of states pi depending on p. Moreover, the closest zero-discord states to 
p are given in terms of the corresponding optimal von Neumann measurement. 

The paper is organized as follows. The definitions of the quantum discords and the Bures distance are given in 
Sect. II, together with their main properties. In Sect. Ill we show that the Bures-GQD for a pure state coincides 
with the geometric measure of entanglement and is simply related to the highest Schmidt coefficient. We explain 
this fact by noting that the closest zero-discord states to a pure state are convex combinations of orthogonal product 
states. The link between the minimal Bures distance to the set of zero-discord states and ambiguous quantum state 
discrimination is explained and proved in the following Sect. IV. Finally, we study in Sec. V the case of two qubits. In 
that section the Bures-GQD is calculated explicitly and compared with the other quantum discords, and we determine 
the closest zero-discord states. The last section contains some conclusive remarks and perspectives. 

II. DEFINITIONS OF THE QUANTUM DISCORDS 

In this work we consider a bipartite quantum system AB with Hilbert space % = "Ha ® "Hb , the spaces "Ha and 
T-Lb of the subsystems A and B having arbitrary finite dimensions ua and ns- The states of AB are given by density 
matrices p on "H (i.e., hermitian positive N x N matrices with unit trace, with TV = UAnB)- The reduced states of A 
and B are defined by partial tracing p over the other subsystem. They are denoted by pA — trB{p) and pB — tiAip)- 

Let us first recall the definition of the quantum discord ]2, 3]. The total correlations of the bipartite system in 
the state p is described by the mutual information Ia:b{p) — S{pa) + S{pb) — S{p), where S'(-) stands for the 
von Neumann entropy. The amount JB\Aip) of classical correlations is given by the maximal reduction of entropy 
of the subsystem B after a von Neumann measurement on A. Such a measurement is described by an orthogonal 
family {nf} of projectors acting on Ha- Hence Jb\a{p) = inax|^^}{S'(p_B) — X^i 9i'5'(Ps|i)}; where the maximum 
is over all von Neumann measurements {tt^}, qi = txiir^ ® \ p) is the probability of the measurement outcome i, 
and pB\i — q^^ tryii(7r/^ (S) 1 /)) is the corresponding conditional state of B. The quantum discord is by definition the 
difference 6a{p) — Ia:b{p) ~ Jb\a{p) between the total and classical correlations. It measures the amount of mutual 
information which is not accessible by local measurements on the subsystem A. Note that it is asymmetric under the 
exchange A ^^ B. It can be shown ]32] that 5a{p) > for any p, and SAic^A-cx) = if and only if 

C^A-c\^^q^\a^){ai\®(TB\^, (1) 

i=l 

where {|ai)}"^]^ is an orthonormal basis of "Ha, <^B\i are some (arbitrary) states of B depending on the index i, 
and q-i > are some probabilities, X^i 9* ~ 1- The fact that Sa{cfa-ci) = follows directly from Ia:b{o'A-ci) = 
S (tr Aic A-ci)) — 'l2i9i^i'^B\i) "£ Jb\a{<^A-ci) and from the non-negativity of the quantum discord. For a bipartite 
system in the state <ja-c1i the subsystem A is in one of the orthogonal states ja^) with probability g,, hence A behaves 
as a classical system. For this reason, we will call A-classical states the zero-discord states of the form (1). We 
denote by Ca the set of all ^-classical states. By using the spectral decompositions of the aB\ii any ^-classical state 
(JA-c\ G Ca can be decomposed as 

CTA-ci = ^^gy|ai)(aj| ® \fij0{Pj\z\ (2) 

where, for any fixed i, {|/3j|i)}?Ji is an orthonormal basis of "Hs, and qij > 0, ^^ ■ qij = 1 (note that the |/?j|i) need 
not be orthogonal for distinct i's). One defines similarly the set Cb of S-classical states, which are the states with 
zero quantum discord when the subsystem B is measured. A state which is both A- and _B-classical possesses an 
eigenbasis {\ai) <Si |/3j)}"^]^"f,]^ of product vectors. It is fully classical, in the sense that a quantum system in this state 
can be "simulated" by a classical apparatus being in the state {i,j) with probability qij. 

Let us point out that Ca, Cb, and the set of classical states C are not convex. Their convex hull is the set S of 
separable states. A state cFscp is separable if it admits a convex decomposition crgcp — '12m 9m|</')n)(0m| (Xi \'4'm){''l'm\, 
where {|(/'m)} and {|V'm)} are (not necessarily orthogonal) families of unit vectors in Ha and Hb and q^ > 0, 
12m1"i ~ 1- ^'^^ pure states, ^-classical, B-classical, classical and separable states all coincide. Actually, according 
to (2) the pure ^-classical (and, similarly, the pure i?-classical) states are product states. 

The geometric quantum discord (GQD) of a state p of AB has been defined in [21] as the square distance of p to 
the set Ca of A-classical states, i.e., 

DA{p)=d2{p,CAY^ min d2{p,aA-c\Y (3) 

"■A-clGC/i 



where d2{p, u) — (tr[((0— cr)^])^/^ is the Hilbert-Schmidt distance. Instead of taking this distance, we use in this article 
the Bures distance 



dB{p,a) = 2{l-^Fi^)) (4) 

where p and a are two density matrices and F{p, a) is their fidehty [1, 14, 33], 

F{p,a) ^ WVpV^WI = [tr([V^pV^i/2)]2 _ (5) 

It is known that (4) defines a Riemannian distance on the convex cone £ C Mat(C,N x N) of all states of AB. Its 
metric coincides (apart from a numerical factor) with the quantum Fisher information which plays an important role 
in quantum metrology [17]. Moreover, ds satisfies the following properties ]1, 33]: for any p, cr, ci, and (T2 G £, 
(i) convexity: if ?7i,??2 > and 771 + 772 = 1, then ^5(^,771(71 +?72CT2) < '7irfB(/5, ^i) + r]2d{p,cr2); 
(ii) ds is preserved by unitary conjugations: if [/ is a unitary operator on T-L, then dsiU pU\ U(tW) = dsip, c); 
(iii) ds is monotonous under the action of completely positive trace preserving maps T from Mat(C, N x N) into itself: 
for any such T, dsiTp^Ta) < dB{p,cr)- 

Note that the Hilbert-Schmidt distance d2 also satisfies (ii) but fails to satisfy (iii). The monotonous Riemannian 
distances on £ have been classified by Petz [16]. 

We define the geometric quantum discord (GQD) as 

Da{p) - dB{p,CA) = 2(1 - VFa{p)) , Fa{p) = max F{p,<jA-ci) ■ (6) 

The property (ii) above implies that Da and Dj^ are invariant under conjugations by local unitaries, p h- > Ua ^ 

UbpUa ® U^, since such transformations leave Ca invariant. By property (iii). Da is monotonous under local 
operations involving von Neumann measurements on A, generalized measurements on B, and classical communication. 
By analogy with (6) one can define two other geometrical measures of quantum correlations: the distance D{p) ~ 
dB{p,C) to the set of classical states C and the geometric measure of entanglement E{p) ~ dB{p,S) = 2(1— ^/Fs{p)), 
where Fg(p) is the maximal fidelity between p and a separable state Csop G S. The second measure satisfies all criteria 
of an entanglement measure ]13] (in particular, it is monotonous under local operations and classical communication 
by the property (iii)). It is related to other entanglement measures defined via a convex roof construction thanks to 
the identity ]18] 

Fs{p)= max Vp^F^d*™)) , (7) 

{p-.}.{|*™>}tr 

where the maximum is taken over all pure state decompositions p — X]mP™l^'")(^™l *^^ P (with ||^m|| = 1 and 
Pm > 0, J2Pm — 1) ^^'^ Fs{\'^„i)) — maxo-e5 -F(|\l/m)(^E'm|,o-). The measure E{p) is a geometrical analog of the 
entanglement of formation ]34]. The latter is defined via a convex roof construction from the von Neumann entropy 
of the reduced state, Eeof{p) = min{p^}_{|^^)} X]m-P™-^EoF(|*m)) with E'Eopd^m)) = S{trB |*m)(*m|). 
Since C C C^ C 5, the three distances are ordered as 

E{p) < Da{p) < D{p) . (8) 

This ordering of quantum correlations is a nice feature of the geometric measures. In contrast, the entanglement of 
formation Eeof{p) can be larger or smaller than the quantum discord Sa{p) [35, 36]. 

III. THE BURES GEOMETRIC QUANTUM DISCORD OF PURE STATES 

We first restrict our attention to pure states, for which one can obtain a simple formula for Da in terms of the 
Schmidt coefficients pi. We recall that any pure state |^) G Ha 'S>'Hb admits a Schmidt decomposition 

n 

I*> = "EVmII^«)®Ix*) (9) 



1=1 

where n — minjriyi, 773} and {|<y5i)}"^i (respectively {|Xj)}?Ji) is an orthonormal basis oH-La {'Hb)- The decomposi- 
tion (9) is unique, the pi and \ipi) (respectively \xj)) being the eigenvalues and eigenvectors of the reduced state pA 
of A (respectively pB of B). 



Theorem 1. //pip = |^)(^| is a pure state, then 

Da{p^) = D{p^)=E{p^) = 2{1-^JI^^) (10) 

where /imax is the largest Schmidt coefficient pi . If moreover a single coefficient pi is equal to /imax then the closest 
A-classical (respectively classical, separable) state to p^, is the pure product state a = \ipma.x (Xi Xmax){^max 'Si Xmaxl? 
where |(/?niax) o-i^d |xmax) «'~e the vectors corresponding to /imax ii^ the decomposition (9). If r > 1 Schmidt coefficients 
Pi are equal to Pmax, say /Umax — Pi = ■ ■ ■ = Pr > Pr+ii ■ ■ ■ , Pn, then infinitely many A-classical (respectively classical, 
separable) states a minimize the distance dB{p^i,<j). These closest states a are convex combinations of the orthogonal 
product states \ai <Si f3i){ai <Si f3i\, I = I, . ■ . ,r, with \ai) = X^Li uul^i) and \/3i) = X][=i u*i\xi) , where iuii)li^-^ is an 
arbitrary r x r unitary matrix and \ipi) and \xi) are the vectors in the decomposition (9). 

The expression (10) of the geometric measure of entanglement E{p^,) is basicaUy known in the hterature, and 
the closest separable states to pure and mixed states have been investigated in [18]. By inspection of (9) and 
(10), Da{p^!) = if and only if |^) is a product state, in agreement with the fact that A-classical pure states are 
product states (the same holds for the other quantum correlation measures D and E). Moreover, from the inequality 
/imax ^ 1/" (following from X]i=i A*« ~ 1) ^''^^ deduces that DA{p^i) < 2(1 — l/^/n). The maximal value of Da is 
reached when pi — \/n for any i, that is, for the maximally entangled states (recall that such states are the pure 
states with maximal entropy S{pa) — ~ X]r=i Mi ^^Pi = ln(n)). Remarkably, the maximally entangled states are the 
pure states having the largest set of closest separable states (this set is a (n^ + n — 1) real-parameter submanifold of 
£). For instance, in the case of two qubits (i.e., for ua = tib = 2), the Bell states |$^) = (|00) ± |ll))/\/2 admit as 
closest separable states the classical states 

a±= ^ g,|az)(an®|A)(/3d , \ai) ^ uoi\Q) + uu\l) , |A) = "w|0) ± <;|1) (H) 

1=0,1 

with UQiUQm + u^jUim — 5ml and (// > 0, go + 91 = 1- Interestingly, typical decoherence processes (such as pure 
phase dephasing) transform p$± into one of its closest separable state (|00)(00| + |ll)(ll|)/2 at times t ^ idee, where 
tdec is the decoherence time. Slower relaxation processes modifying the populations in the states |00) and |11) do 
not further increase the distance to the initial state p$± . The situation is different for a partially entangled state 
1^) = \/7^|00) + a/mTIH) with pi > pq: then the closest separable state is the pure state |11), but |^) evolves 
asymptotically to a statistical mixture of |00) and |11) when the qubits are coupled e.g. to thermal baths at positive 
temperatures. 

Proof. For a pure state p.p, the fidelity reads F{p<ii,a-A-c\) = (^Ica-ciI^)- Replacing o-a-cI in (6) by the right-hand 
side of (2) we get 

Fa{p^)= ^^max {^ g,,|(«, ® /3,|,|vl/)|2| ^ max |(« ® /3|*)P (12) 

where we have used "^^^qij ~ 1. Thanks to (9) and to the Cauchy-Schwarz inequality, for any normalized vectors 
|a) e Ha and |/3) € He one has 

|(a0/3|*)| <^7^|(a|<^,)(/3|x,)| < VAW;^|(a|<^^)(/3|x^)| 

i i 

i i 

Let us first assume that pi — /imax > P2, • • • , Pn- Then \{a <Si /3|^)| = ^//^max if and only if |a) ~ \ipi) and |/3) ~ \xi)- 
Thus the maximal fidelity between p^j, and an A-classical state is simply given by the largest Schmidt coefficient /imax 
and (10) follows from (6). Moreover, the maximum in the second member of Eq.(12) is reached when a single qij is 
non-vanishing, say qij — SuSji, and |ai) — \(pi), |/?i|i) — \xi)- This means that the closest A-classical state to pq, is 
the pure product state \(pi (g) Xi){'Pi ^ Xi\- 

More generally, let pi = ... = pr = /imax > Pr+i, ■ ■ ■ , Pn- In order that the first and second inequalities in 
the first line of (13) be equalities, arg((Q;|(^i)(/3|xi)) = must be independent of i and either |a) belongs to Imax — 
span{|</?i), . . . , \(pr)} or |/3) belongs to span{|xi), . . . , |xr)}- Furthermore, the Cauchy-Schwarz inequality in the second 
line is an equality if and only if |(a|(/?i)| ~ \{P\Xi)\ for all i. This shows that from any orthonormal family {|q!;)}[^]^ 
of V^ax, one can construct r orthogonal vectors \ai (E) /?/) satisfying \{ai (g) /3/|^)| = ^/imax for alH = 1, . . . ,r, with 
il^lxi) = 6 ((y9;|a). The probabilities {qij} maximizing the sum inside the brackets in (12) are given by qij — qi if 



i ^ i "£ r and zero otherwise, where {qiYi=i is an arbitrary set of probabilities. The corresponding A-classical states 
with maximal fidelities F{pq,,a) are the classical states a = X]i=i lil'^i ® Pi){'^i ® A|- 

The equality between the correlation measures Da^ D, and _E is a consequence of the fact that the closest states to 
pq, are classical states. Such an equality is reminiscent from the equality between the entanglement of formation -EeoF 
and the quantum discord Sa for pure states. Let us notice that it does not hold for the Hilbert-Schmidt distance, for 
which the closest A-classical state to a pure state is in general a mixed state. Actually, one infers from the expression 

d2ip^,<yA-cif - tr[(|*>(*| - c7A-c\f] = 1 - 2F(p*,aA-ci) +tr((7l,i) (14) 

that the closest A-classical state results from a competition between the maximization of the fidelity F{p^, cr^-ci) and 
the minimization of the trace tr((T^_^j), which is maximum for pure states. For instance, one can show that the closest 
A-classical state to the Bell states |$^) are mixed two-qubit states, see Sec VD. The validity of theorem 1 is one of 
the major advantage of the Bures-GQD over the Hilbert-Schmidt-GQD. 

IV. THE BURES GEOMETRIC QUANTUM DISCORD OF MIXED STATES 

A. Link with minimal error quantum state discrimination 

The determination of Da{p) is much more involved for mixed states p than for pure states. We show in this 
section that this problem is related to ambiguous quantum state discrimination (QSD). As it has been recalled in 
the introduction, in ambiguous QSD a state pi drawn from a known family {pi}"^^ 'with prior probabilities {r/i}"^^ ^s 
sent to a receiver. The task of the latter is to determine which state he has been given with a maximal probability of 
success. To do so, he performs a generalized measurement and concludes that the state is pj when his measurement 
result is j. The generalized measurement is given by a family of positive operators Mi > satisfying "^^Mi ~ 1 
(POVM). The probability to find the result j is Pj\i ^ ti{Mjpi) if the system is in the state pi. The maximal success 
probability of the receiver reads 

UA 

Ps'^HlPi^Vt}) = niax V?7,tr(M,p,) . (15) 

POVM {Mi} '^ — ' 

2 — 1 

Theorem 2. Let p be a state of the bipartite system AB with Hilbert space % — "Ha ® "Hb and a = {\ai)}^^i be an 
orthonormal basis o/Ha- Consider the subset of A-classical states with a diagonal marginal in the basis a, 

CA{ct) = {o'A-ci E Ca T ct is an eigenbasis of trBif^A-ci)} ■ (16) 

Then the maximal fidelity F{p,CA{ct)) — max F{p,aA-ci) of p to this subset is equal to 

■n.A 

F{p,CA{a)) = Ps^'^■''■{{p^,Vi\) = max^r;,tr(n,p,) (17) 

*- '-" i=i 

where Pg^ ^' '{{pi,rii}) is the maximal success probability over all von Neumann measurements given by orthogonal 
projectors Hi of rank ns (that is, self-adjoint operators on % satisfying HiHj = SijHi and dim(ni'H) = ns), and 

Vi ^ {ai\pA\ai) , Pi = 'n7^\/p\a-i){ai\®ly/p (18) 

(if rji = then pi is not defined but does not contribute to the sum in (17)). 

This theorem will be proven in Sec. IV B. Note that the pi are quantum states of AB if rji > 0, because the 
right-hand side of the last identity in (18) is a nonnegative operator and rji is chosen such that tr(pi) ~ 1. Moreover, 
{Vi}7=i is a set of probabilities (since rji > and X^i '?« = ^''^(p) ~ 1) and {pi,??i}"=i defines a convex decomposition 
of P, P^Y^iViPi- 

Let us assume that p is invertible. Then the application of a result by Eldar [37] shows that the POVM maximizing 
the success probability Ps{{pi,rii}) in (15) is a von Neumann measurement with projectors Hi of rank ub, i.e., 

F{p,CA{a))^P°^'-■''■{{p„V^})^Ps°''\{p^,V^}) , P>0. (19) 

In fact, one may first notice that all matrices pi have rank r^ = ub (in fact, pi has the same rank as riip~^''^pi = 
\oii){ai\ ® l-y/p and the latter matrix has rank ns)- Next, we argue that the pi are linearly independent, in the sense 



that their eigenvectors \£,ij) form a Hnearly independent family {|Cji)}i=i' "^ of vectors in H. Actually, a necessary 
and sufficient condition for \^ij) to be an eigenvector of pi with eigenvalue Ay > is j^y) — {Xijrii)^^y/p\ai) ® \Cij), 
\Qj) € Hb being an eigenvector of Ri ~ {ai\p\ai) with eigenvalue Xijrji > 0. For any i, the hermitian invertible matrix 
Ri admits an orthonormal eigenbasis {\Cij)}^=i- Thanks to the invertibility of y/p, {\£,ij)}lZi 'n'^ is a basis of H and 
thus the states pi are linearly independent. It is shown in [37] that for such a family of linearly independent states 
the second equality in (19) holds true. 

The following result on the Bures-GQD of mixed states is a direct consequence of Theorem 2. 

Theorem 3. For any state p of the bipartite quantum system AB, the fidelity to the closest A-classical state is given 
by 

n,A 

Fa{p) = max maxy^ tr[ni^|ai)(aj| lJp\ (20) 

{la,}} {n,} ^^-^ L Vi /x I vrj 

«— 1 

where the maxima are over all orthonormal basis {\oii)} of "Ha o,nd all families of orthogonal projectors {Hi}"^]^ of 
Ha $5 T~{-B with rank ub ■ Hence, using the notation of Theorem 2, 

F^(p)=maxP7'--^-({p„,7j). (21) 

If p > then one can replace Pg^ ^' ' in (21) by the maximal success probability (15) over all POVMs. 

It is noteworthy to observe that the basis vectors \ai) can be recovered from the states pi and probabilities rji by 
forming the square root measurement operators Mi = rjip~^''^pip~^''^, with p = 'Yli'HiPi {^^ assume here p > 0). 
Actually, such measurement operators are equal to the rank-n^ projectors Mi ~ \ai){ai\ ® 1. By bounding from 
below Pg^ ^' '({/9i,?7i}) by the success probability corresponding to IIj = Mi, we obtain 

"A 

Fa{p) > max V trs [(a»|Vp|a,)'] • (22) 

{!".)} -^ 

The square root measurement plays an important role in the discrimination of almost orthogonal states [38] and of 
ensembles of states with certain symmetries [27[. 

To illustrate our result, let us study the ambiguous QSD task for some specific states p. 

(i) If p is an A-classical state, i.e., if it admits the decomposition (1), then the basis {|ai)} maximizing the optimal 
success probability in (21) coincides with the basis appearing in this decomposition. With this choice, one obtains 
r]i = Qi and pi = \ai){ai\ <Si (7B\i for all i such that qi > 0. The states pi are orthogonal and can thus be perfectly 
discriminated by von Neumann measurements, so that Fa{p) ~ -Pj^ ^' '{{Pi,Vi}) = 1- Reciprocally, if Fa{p) = 1 
then Pg^ ^' (p) = 1 for some basis {[oi)} of Ha and the corresponding pi must be orthogonal, that is, pi ~ liipilii 
for some family {11^} of orthogonal projectors with rank ub- Hence p — ^^ r]iPi — ^^ "qiHipiHi, ^ = ^^ Iliy^IIi, 
and (18) entails r\iPi = r\i^iPi^i = ^\oL^{ai\ ®\^ = ^^i\oLi){ai\® 111^^, implying 11^ = \ai]{pii\® 1. Thus p is 
A-classical. This was of course to be expected since Da{p) = if and only if p is A-classical, see Sec. II. Therefore, we 
can interpret our result (21) as follows: the non-zero discord states p are such that the states (18) are non- orthogonal 
and thus cannot be perfectly discriminated for any orthonormal basis {\ai)} ofHA- 

(ii) li p ~ p^ij is a pure state, then all pi with r/i > are identical and equal to pij,, so that Pg^ = Pg^ ^' ' = 
suprjj.i ^j ?7i(^|ni|^) = r/max- One gets back the result Fa{p^) ~ Mmax of Sec. Ill by optimization over the basis 

(iii) Let us determine the states p having the highest possible GQD, i.e., the smallest possible fidelity Fa{p)- 

Proposition. If ha < ns, the smallest fidelity Fa{p) for all states p of AB is equal to I/ua- IfruA < ns < {r + l)nA 
with r — 1,2, . . ., the states p with Fa{p) ~ ^/nA are convex combinations of the r maximally entangled pure states 

l*fe> - n-/^' E7A l#^) ® IV'f ^), k^l,...,r, wtth (#Vf ) = S.J and (^f Vf > = 5kA,. 

We deduce from this result that the GQD Da{p) varies between and 2 — 2/y/riA if n-A < "-b- By virtue of 
Theorem 1, the proposition above, and the inequality E{p) < Da{p)) the geometric measure of entanglement E{p) 
also varies between these two values. This means that the most distant states from the set of A-classical states Ca are 
also the most distant from the set of separable states S. If ua < nB < 'in a, these most distant states are maximally 
entangled pure states. 



Proof. The success probability Pg^ ^' ' must be clearly larger than the highest prior probability rj^ax — inaxi{?7i}. 
(A receiver would obtain Pg = r^max by simply guessing that his state is Pi^^^, with 77^^,^ = r/max; a. better strategy 
is of course to perform the von Neumann measurement {11^} such that H^^^^ projects on a ris-dimensional subspace 
containing the range of pi^^^ ; this range has a dimension < ns by a similar argument as in the discussion following 
Theorem 2). In view of (21) and by using ?7max > ^/"^a (since ^^ rji = 1) we get 

Fa{p) > — (23) 

riA 

for any mixed state p. 

When HA < ns the bound (23) is optimum, the value I/ua being reached for the maximally entangled pure states 
p = p^, see Sec. III. Thus I/ua is the smallest possible fidelity. Let p be a state having such a fidelity Fa{p) = l/nA- 
According to (21) and since it has been argued before that Pg^ ^' ' > ?7inax ^ ^/nA, Fa{p) = l/nA means that 
pop V. ({p^^Ty^}) — I/ua whatever the orthonormal basis {|ai)}. It is intuitively clear that this can happen only if 
the receiver gets a collection of identical states pi with equal prior probabilities r]i = I/ua- A rigorous proof of this fact 
is given in the Appendix. From (18) and p = "^rjipi we then obtain {ai\pA\Q:i) = I/ua and Pi — p for any i = 1, . . . , ua 
and any orthonormal basis {|ai)}. The first equality implies pA ~ ^/tia- By replacing the spectral decomposition 
p = Z]fePfe|*fe)(*fc| of P into (18), the second equality yields tr_B(|*fe)(\I';|) = n'^Ski for aU fc, I with p^p/ 7^ 0. Taking 
advantage of this identity for k ~ I, one finds that the eigenvectors I'f/j,) of p with positive eigenvalues p^. have all 

their Schmidt coefficients equal to I/ua, i.e., their Schmidt decompositions read j^^,) — rij^ ' 'J27=i I'^i ) '^ l^i )■ 
Moreover, trB(|*fe)(*,|) = is equivalent to V^j'^'-LVJ^' with v{j''' = span{|V'f ^); i = l,...,n^} C Hb- If ng < 
(r + l)nA then at most r subspaces Vg may be pairwise orthogonal. Thus at most r eigenvalues pk are nonzero. 

Let us briefiy discuss the case ua > ng. In that case the smallest value of the maximal fidelity Fg{p) to a separable 
state (see Sec. II) is equal to l/ng and Fg{p) — I/ub if and only if p is a pure product state. This is a consequence 
of (7) and of the bound Fs{p<i,) > l/ns for pure states p^. As a result, E{p) varies between and 2 — 2/y/n with 
n = inm{nA,nB}- We could not establish an analog result for the GQD Da{p)- The following reasoning shows that 
the bound (23) is not optimal, i.e., there are no states p with fidelities Fa{p) equal or close to I/ua when ua > ns- 
Indeed, by generalizing the proof above, such states should be close to pure states. Since for ua > nB pure states 
have fidelities larger than I/ub (see Sect. Ill), one may expect that a state close to a pure state has a fidelity close 
to I/hb or larger. This can be shown rigorously by invoking the inequality Fa{p) > ||p||/nB, where ||p|| = max{pi:} 
is the norm of p. This inequality is obtained by bounding Fa{p) by the maximal fidelity to the set of pure product 
states, 

Fa{p) > ,, ,max (« (g> p\p\a (g> l3) > max{pkp^^lj > !^^^^^M , (24) 

||a|| = ||/3||=l k UB 

where we have used in the second inequality the argument in the proof of Theorem 1 to express max||Q,||^||^||^i |(a <S) 
/3|vl'fc)p as the largest Schmidt coefficient /Xmax of |^fc). Note that all bounds in (24) are optimal for pure states only. 

B. Derivation of the variational formula (20) 

To prove Theorems 2 and 3, we start by evaluating the trace norm in (5) by means of the formula ||r||i = 
max[/ 1 tr([/X)|, the maximum being over all unitary operators on "H. Using also (2), one gets 



\/F(p, <TA-c\ = max|tr(C/^^crA-ci) | 

In the third line we have replaced the maximum over unitaries [/ by a maximum over all orthonormal basis {|$ij)} 
of "H (with |$ij) = U^ai (X) /3j|i))- The last equality in (25) can be explained as follows. The expression in the last 



u 



max I 

u 



line is clearly larger than that of the third line; since for any i and j one can choose the phase factors of the vectors 
|$y) in such a way that {^ij\y/p\ai (E) f3j\i) > 0, the two expressions are in fact equal. 

One has to maximize the last member of (25) over all orthonormal basis {|/3j|i)} of Hb and all set of probabilities 
{q-ij}- The maximum over the probabilities qij is easy to evaluate by using the Cauchy-Schwarz inequality and 
X^i lij = 1- It is reached for 

We thus obtain 

F{p,CAia)) = max max F{p, a a-ci) = max max V|(*j|j|/3j|j)| (27) 

{\P3\^)} ili,} {|/3j|i)} {|*ij)}77 

where we have set |V'j|i) = {o^il^/pl^ij) ^ ^s- We proceed to evaluate the maximum over {|/3j|i)} and {|$ij)}- Let 
us fix i and consider the orthogonal family of projectors of "H of rank ub defined by 

n, = Y,\<i>,,){^,,\. (28) 

j 
By the Cauchy-Schwarz inequality, for any fixed «, one has 

rf?«^^,IlK^j|*l'^j|*)r^IIll^j|*ll'=*^[n.Vp|a.)(a.|®lVp] . (29) 

Note that (29) is an inequality if the vectors iV'jii) are orthogonal for different j's. We now show that this is the 
case provided that the |$ij) are chosen appropriately. In fact, let us take an arbitrary orthonormal basis {|$ij)} of 

7i and consider the hermitian ub x ub matrix S^^' with coefficients given by the scalar products S*-^ — {'4'j\i\'4'k\i)- 
One can find a unitary matrix F*^*^ such that S'^^^ — (y(*))t5'(0]/(0 ig diagonal and has nonzero diagonal elements 
in the first r^ raws, where r^ is the rank of S'^^K Let |$ij) = X]I=fi ^u l^u)- Then {|^ij)} is an orthonormal basis 
of H and ^ - |$ij)($ij| — Hi. Moreover, the vectors \'ipj\i) — (Q:i|y/p|$ij) = X)j*fi 1^1 IV';|i) form an orthogonal set 
{\^j\i)}''i=i ^^'^ vanish for j > r^. Therefore, for any fixed orthogonal family {Hj}"^]^ of projectors of rank ub, there 
exists an orthonormal basis {|$ij)} of H such that (28) holds and the inequality in (29) is an equality. Substituting 
this equality into (27), one finds 

F{p,CA{a)) = max^ ll^j|i|P = maxY^ tr[n,^|Q;,)(o;i| ® 1 ^] , (30) 

i,j ' i 

which yields the result (17). The formula (20) is obtained by maximization over the basis {jofi)}. 

C. Closest ^-classical states 

The proof of the previous subsection also gives an algorithm to find the closest A-classical state to a given mixed 
state p. To this end, one must find the orthonormal basis {|q;°^ )}r=i "-"^ '^a maximizing Pg'^ ^' '({pi,77i}) m (21) 
and the optimal von Neumann measurement {11°^ }7=i yielding the minimal error in the discrimination of the states 
p°^ with probabilities 77°^ which are associated to {\a°^ )} by the formula (18). 

Theorem 4. The closest A-classical state to p is 

'^p = j^) E i"r)(«r I ® (ar iv^nr vpi^Da (31) 

where {\a°/')}'l^, and {n°P'}^4 are such that Fa{p) = l:^^AK^' ^\aT){^T\ ® 1 Vp\' «ee Eq.(20). 

Proof. This proof follows directly from that of the previous subsection. Actually, by using the expression (26) 
of the optimal probabilities qij and the fact that the Cauchy-Schwarz inequality (29) is an equality if and only if 
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\Pj\i) — IV'j|i)/l!V'i|i|l when HV'jiill ¥" Oj we conclude that the closest A-classical state to p is given by Eq.(2) with 

^ E"filV'j|i)(V'j|i| 

9«aB|« - 2^9yl/3jK)(/?jKl = ~ (32) 

The denominator is equal to Fa{p), see (30). The numerator is the same as the second factor in the right-hand side 
of (31) as by construction IV-jIj) = {a°/^\y/p\^ij) and J2^ |$y)($ij| = 11°^*. 

V. THE TWO-QUBIT CASE 

A. Maximal success probability 

If the subsystem A is a qubit, the set {pi} contains only ua = 2 states. The minimal error and optimal measurements 
to discriminate the states pi with prior probabilities rji are then easy to determine [24] . In our case we must maximize 
the success probability over all von Neumann measurements given by orthogonal projectors IIq and Hi having the 
same rank ns- One starts by writing the projector Hi as 1 — Ho in the expression of the success probability, 

P^'H{P^,V^}) = ^0 tr(nopo) + VI tr((l - no)pi) = i(l - tr A) + tr(noA) (33) 

with A = rjopo — ryipi. Thanks to the min-max theorem, the maximum of tr(noA) over all projectors Hq of rank ns is 
equal to the sum of the ub highest eigenvalues Ai > • • • > A„g of the hermitian matrix A, and the optimal projector 
IIq^ is the spectral projector associated to these highest eigenvalues. One has 

1 "B 

^s°'*'-''-({P«-^J) =nmxPf '>({p„ryj) = -(l-trA) +^A, . (34) 

^ '-'" 1=1 

For the states pi and probabilities rji defined by (18), A — ^/p{\aQ){aQ\ — |ai)(ai|) (g) 1 y/p where {|ai)}i=o i^ ^^ 
orthonormal basis of C^. For any such basis the operator inside the parenthesis in the last formula is equal to 
(Tu = X]m=i ^nifm for somc Unit vector u € R"^, where cti, (T2, and a^ are the Pauli matrices. Reciprocally, one can 
associate to any unit vector u e R'^ the eigenbasis {|ai)}i^o ^^ '''"• Hence, by substituting (34) into (21) we obtain 

1 "-B 

F^(p) = - max|l-trA(u) + 2VA;(u)| , A{u) ^ y^a^® 1 ^^ , (35) 

Z Hull — 1 L ^ — ^ J 

where A;(u) are the eigenvalues of A(u) in decreasing order. The term 1 — trA(u) is equal to tr(p^(l — a^)) = 
2(ai|p^|Q;i). Note that —p < A(u) < p, so that X]/!fi ^li'^) ^ YM=iPh where p; are the eigenvalues of p in decreasing 
order. If p > then A(u) has ub positive and ub negative eigenvalues; in such a case (34) reduces to the well known 
expression [24] P^'^-^dft, ,7,}) = P7'({ft,,?J) -(1 + tr|A(u)|)/2. 

In order to evaluate Fa{p) one needs to diagonalize the 2nB x 2nB matrix A(u) and to determine the maximum 
over the unit vectors u in (35). We do this in the next subsection for a special class of two-qubit states. 

B. Bures-GQD for states with maximally disordered marginals 

With the aim of comparing Da with the other quantum discords, let us consider a simple family of two-qubit states 
for which all discords can be easily calculated: the states p with maximally disordered marginals pa = pB = 1/2. 
Any such state can be written up to a conjugation by a local unitary Ua ® Ub (which leaves all discords unchanged) 
as 139] 

1 ^ 

/5= -fl® 1 + ^ c„cr„ (8)cr„j (36) 

r?i— 1 

where the real vector c = (ci, C2, C3) belongs to the tetrahedron with vertices F± = (±1, =f1, 1) and G± = (±1, ±1, —1), 
that is, 

Po = t(1 -ci - C2 - C3) > , p„ = -(1 + C1+C2 + C3- 2cm) > , m = 1,2,3 . (37) 
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We now calculate Da{p) for two qubits in the state (36). In the standard basis {|00), |10), |01), |11)}, the matrices 
p and y/p are given by 



/ 1 + C3 Cl - C2 \ 

I-C3 C1+C2 

Ci + C2 1 - C3 

V Ci - C2 1 + C3 / 



v~p 



( to + h ti~t2\ 

to - h ti+h 

h+t2 to-t3 

\t1-t2 to + 13/ 



(38) 



where 



4to = \/P0 + VPI + \/P2 + \/P3 , ^ti = -y/Po- y/Pl + VP^+ VP3 , 



(39) 



with similar formulas for t2 and t^ obtained by permutation of the indices 1,2,3. The real parameters ti, satisfy 
4(t§ + ^^+^2+^3) = ti p = 1. Let 9 and (p be the angles defined by the vector u = (sin 6 cos (f), sin 6 sin (p, cos 6) on the 
unit sphere. In the standard basis the matrix A(u) = ^/pcTu ® l-^/p reads 



A(u) = i 



f{a3 + b3)cose Qsin0 C^sin^ 

C^sin^ (03 — 63) cos 

^0sin(^ (— a3 + 63)cos0 





Qsine 
C;sin0 



\ 



V 



, sin(^ 



(40) 



{-as - 63) cos 0/ 



^0 sin 9 

with ^0 — fli cos (/i + ia2 sin 0, ^^ = 61 cos + 162 sin 0, and 

&m = 8{tl + 4J - 1 , ai = 8(ioti + t2t3) , 02 = 8{tot2 + hts) , 03 = 8(toi3 + ^1*2) 
One finds 

as = 2(-VPoP3 + V-P1P2) , &3 ==2(^poP3 + \/-PlP2) • 

Similar expressions holds for the other coefficients ai, 02, &i, and 62 by permuting the indices 1, 2, 3. 
The eigenvalues of A(u) come in opposite pairs (A±(u), — A±(u)) with 



1 



A±(u) = ^ ^Jblcos'^9+\C^\'^smH±^alcos'^9+\^4,\'^sm'^0 



(41) 
(42) 

(43) 



in agreement with tr A(u) = 1 — 2{ai\pA\oii) = 0. But \am\ < &m by (42), thus the second square root is smaller than 
the first one. One deduces from (35) that 



Fa{p) 



1 + max -i/6§ cos2 9 + {bf cos^ (j) + bl sin 0) sin 9 

2 \ 9.(p 



(44) 



If the b„i are distinct from each other, the maximum over the angles 9 and (p is reached for cos^ 9, cos^ (p E {0, 1}, that 
is, for u = ±63;, ±ej,, or ie^. If 61 = &2 > ^3, this maximum is reached for cos^ 9 = and cose/) arbitrary. Thus 



Fa{p) = 



- and Da{p) = 2(i- ^ 



(45) 



with 



max {6m} = - max| ^^(1 + ci)^ - (c2 - 03)^ + \/{l - c\f - [02 + 03)^ , 
1=1,. ..,3 Z I 



V/(1 + C2)2 - (Ci - C3)2 + v/(l - C2)2 - (ci + 03)2 , 
V(l + C3)2 - (Ci - C2)2 + v/(l - C3)2 - {c, + 02^} 



(46) 



We notice that if \cm\ is maximum for m = mmax then this is also true for &,„, i.e., the maximal number inside the 
brackets in (46) is the (TOmax)th member. Actually, one can show from (42) that 



-"k ' 



m, fc = 1, . . . , 3 



(47) 



12 



(a) 




J C2 



(b) 




K 






/C 


^\ 




/ 


/ \'T 


/<^,^ 


\ '^2 


l\ 




I 


^J 


(c) 




'^l 





l-'-2 




FIG. 1. (Color online) (a) Tetrahedron T with vertices F± = (±1, =1=1,1) and G± = (±1,±1,— 1). The vectors c inside T 
represent physical states p(c) with maximally disordered marginals. The A-classical states in T are located on the segments 
\IK\ \JL\ and [MiV] (thick red dashed lines), with 1(1,0,0), J(0, 1,0), A'(-1,0,0), L(0, -1,0), M(0,0, 1), and iV(0,0, -1). 
The shaded region delimited by the octahedron IJKLMN corresponds to separable states, (b) Cut view of T in the plane 
C3 = 0. The state p(j, ^,0) has a unique closest ^-classical state crp(0, S2,0) with S2 = (4 + -\/35 — \/3)/(8 + \/35 + \/3) — 0.523 
(see (57)). States on the broken lines are at equal Bures distance from the segments [IK] and [JL]. They have infinitely 
many closest A-classical states given by (60). On these lines Da{p) is not differentiable. We show the two closest states of 
the form (36) to the state r(— i,— i,0) (they are obtained by taking </!> = and 7r/2 in (60)): 0-^(81,0,0) and (Tt-(0, si,0) with 
si = — (2 + \/6 — \/2)/(2 + \/Q-\- \/2) ~ 0.518. (c) Cut view of T in the plane Ci = C2. States in the shaded region are separable. 
The state p(|, ^, — |) has a unique closest A-classical state crp(0, 0, S3), where S3 = — (5 + \/33)/(9 + \/33) ~ —0.729. The states 
inside the triangles H^OG± (here 0(0, 0, 0) is the origin) admit infinitely many closest ^-classical states given by (60). States 
belonging to the blue square and triangle in panels (b) and (c) and to a face of T in panel (a) have infinitely many closest 
states given by Eqs.(54), (56), and (55). 



It is clear on (46) that the vectors c such that b takes the highest possible value 6=1 have a single non-vanishing 
component c^. Therefore, in agreement with Ref. [21], the A-classical states with maximally disordered marginals 
are situated (up to local unitary equivalence) on the three segments of the coordinate axes inside the tetrahedron, see 
Fig. 1(a). As a consequence, the separable states with maximally disordered marginals are located inside an octahedron 
with vertices on the three coordinate axes at a distance unity from the origin. This result is well-known [39]. It follows 
here immediately from the fact that such separable states are convex combinations of ^-classical states. For indeed, 
the octahedron is the convex hull of the three aforementioned segments. 

The states (36) with the highest GQD I?^ = 2 — a/2 have 61 = 62 = &3 = 0. In agreement with the proposition of 
Sec. IV A, these states are the four Bell states |$±) = (|00) ± |11))/V2 and |*±) = (|01) ± |10))/V2 located at the 
vertices F± and G± of the tetrahedron (see Fig. 1(a)). Note that an analogous result holds for the Hilbert-Schmidt- 
GQD [21[. 

Let us remark that the inequality (22) is in fact an equality for two-qubit states p with maximally disordered 
marginals pA = PB = 1/2. Actually, for such states one has equal prior probabilities rji = {ai\pA\o:i) ~ 1/2 in the 
QSD problem of Theorem 2. Moreover, it follows from (36) that p is invariant under conjugation by the spin-flip 
operators ai (8) ui and a2 <E) U2- But it has been observed above that an optimal direction u°p* is given by one of the 
coordinate vectors e^ or its opposite. Disregarding irrelevant phase factors, the two eigenvectors jap^ ) and |a°'' ) of 
cru°p' are transformed one into another by the spin- flip operator (7^, with to = 1 if u°p' — ±ej,, ±e^ and/or to = 2 if 
u°P' = ±63:, ie^. Hence the states pi are related by a unitary conjugation, pi = o'm®<^m Pa <^m®'^m- The square root 
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FIG. 2. (Color online) Comparison of the various normalized quantum discords and geometric measure of entanglement for 
a Werner state pw = ^ I + ci*~)(*~| with |*~) = (jOl) - |10))/\/2 (such a state is located on the segment [OG-] in 
Fig. 1(a)). From top to bottom: quantum discord Sa (black boxes), Hilbert-Schmidt-GQD D)^ (red crosses), Bures-GQD Da 
(blue circles), and geometric measure of entanglement E (green diamonds). 

measurement is known to be optimal to discriminate such symmetric set of states with equal prior probabilities [27]. 
Thus Fa{p) is given by the right-hand side of Eq.(22). 



C. Comparison with other quantum discords 



We can now compare Da{p) for two-qubit states p of the form (36) with the other quantum discords defined in 
Sec. II. The corresponding quantum discord 



^a{p) = '^Pu^^2Pu + 2 



\c\ 



ln2(l-|c|) 



■ln2(l + |c|) 



(48) 



i/=0 



has been calculated in [40]. Here, 



maxm I Cm I and the probabilities pi, are given by (37). The Hilbert-Schmidt- 



GQD defined by Eq.(3) is easy to determine for arbitrary two-qubit states [21]; for the particular states (36) it reads 



the GQDs as 



c2 - 



|c| )/4. Since the maximal values of Da and D^ 



(2) 



are 2 — v2 and 1/2, respectively, we renormalize 



Da{p) = 



Da{p) 

2-\/2 



S(2) 



(2)/ 



D'^>{p)=2D'^>{p) 



(49) 



2- V2(l 

1 + C3,|ci -I-C2I 



We also compare Da{p) with the geometric measure of entanglement, given for two qubits by E{p) 
\Jl — C{pYY/^ ]18] where C{p) is the Wootters concurrence. In our case C{p) — 2max{|ci — c; 
1 — C3, 0}. The geometric measure of entanglement is normalized as E{p) = E{p)/{2 — ^/2). 

In Fig. 2, the four quantum correlation measures Da, <5a, Da\ and E are plotted together for the Werner states, 
obtained by taking ci = C2 = C3 = — c with < c < 1 in (36). 

Fig. 3(a) displays the difference Da{p) — 5a{p) for a family of states with maximally disordered marginals depending 
on two parameters. We observe on this figure that 



5a{p) > Da{p) 



(50) 



In contrast, one sees in Figs. 2 and 3(b) that the Hilbert-Schmidt-GQD can be either smaller or larger than 6a, as 
it has been noted previously in other works [41]. We have looked numerically for vectors c inside the tetrahedron 
violating (50) and have not found any such vector. This indicates that (50) is true for any state p with maximally 
disordered marginals. It would be of interest to know if it also holds true for more general two-qubit states. 
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FIG. 3. (Color online) (a) Difference Sa{p) — Da(p) between the quantum discord and the normalized Bures-GQD for states 
with maximally disordered marginals with vectors c inside the triangle of vertices F+, F-, and A'^ (see Fig. 1). The parameters 
u and v,Q<u<l — v<l, are chosen such that ci = — 2u — v -\- 1 = — C2 and C3 = —2v -\- 1. (b) same for the difference 
5a{p) — Dj^ (p) between the quantum discord and the normalized Hilbert-Schmidt-GQD. 



D. Closest ^-classical states 



We now turn to the problem of finding the ^-classical two-qubit states ap satisfying 



dB{p,(Tp) = Da{p) = min ds(p,crA-ci) 



(51) 



These closest A-classical states to p provide useful information on the density matrix p and on the geometry of the 
set of quantum states. Such information is not contained in the discord Da{p) and can be of interest when studying 
dynamical evolutions, for instance, when the two qubits are coupled to their environments and undergo decoherence 
processes (see Sec. III). Like in the two previous subsections, we focus on states p with maximally disordered marginals 
given by Eq.(36), which form a 3-parameter subset T of the (15-parameter) set 8 of all two-qubit states. The geometric 
representation of the states p G T by vectors c in the tetrahedron will give some insight to our results. Moreover, T 
contains many states of interest: maximally entangled states (the four Bell states on the vertices of the tetrahedron), 
the Werner states, separable states (located inside the octahedron), and ^-classical states on the coordinate axes. 
The closest A-classical states for the Hilbert-Schmidt distance d2 are also easy to determine for states p € T by using 
the results of [21]. We will thus be able to compare the geometries on T associated to the two distances ds and d2- 

According to Theorem 4, we must find the optimal basis {\a°^ )} of Ha yielding the maximal value of the optimal 
success probability in (21) and the corresponding optimal measurement operators IIq'' and 11°'' — 1 — IIq^ . Changing 
the order of the vectors in the basis {\a°^ )} amounts to exchange (pq^ , 77q^ ) o (p°^ , 77°'' ) and IIq'' -^ 11°''*. This 
does clearly not modify the optimal success probability and the A-classical state (31). But {|a°'' )} is the eigenbasis 
of the spin matrix ctu for the optimizing direction u — u°p*. Changing the order of the basis vectors thus corresponds 
to inverting u°p'. Therefore, we need only to consider optimal directions u°p* satisfying u°p* • Gm > for m — 1,2, or 
3. Thanks to the results of Sec. VB, we know that if all components Cm of c have distinct moduli, then u°p* — e,,,,^^^, 
where TOmax is the index m associated to the highest modulus |cm|- 

Let us first study the case jcaj > |ci|, |c2|. Then u°p* = 63 and \a°^ ) = \i), {\i)}}^Q being the standard basis of 
Ha — C^. We have seen in Sees. VB and IV A that the states p°^ are optimally discriminated by the square root 
measurement given by the projectors 11°'' = \a°^ ){<^'i^ | <8) 1- In order to find all A-classical states ap satisfying (51) 
we must, however, find all von Neumann measurements {11°'' } maximizing the success probability Ps{{p°^ ,77°'' }). 
As stressed in Sec. V A, Ilg''* is the projector on the direct sum V+ © V_ of the eigenspaces V+ and V_ associated to 
the eigenvalues A+ > A_ > of the matrix A, provided that X± are non-degenerated. If A+ > is simple and A_ = 
is two-fold degenerated, any projector on V+ ® V with V C V_ a one-dimensional subspace of V_ defines an optimal 



projector Ilg'' . By setting 9 = in (40) we find 
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)1 if63>|a3| 

n°P*= <( |00)(00| + |$)($| if63 = a3>0, with |$) € span{|10), |01)} (52) 

|01)(01| + |$')(*'l if h = -as > 0, with 1$') e span{|00), |11)}. 

The condition 63 > jasj is achieved when all p^ are nonzero, see (42). Then the optimal measurement {11°^ } is 
unique and thus p has a unique closest A-classical state. The condition 63 = 03 (respectively 63 — —a^) corresponds 
to P0P3 = (respectively piP2 = 0), i.e., to a vector c belonging to one of the two faces F+F^G± (respectively 
G+G-F±) of the tetrahedron. In this situation one has infinitely many measurements maximizing Ps{{Pi^ ,Vi'^ })• 
Note that (42) implies 63 > |a3|. Moreover, poPs a-nd piP2 cannot both vanish, since otherwise one would have 
bs = 0,3 — 0, in contradiction with our hypothesis |c3| > |ci|, |c2| (which is equivalent to 63 > 61,62 > by (47)). 
Replacing \a°^ ) = \i) and (52) into the expression (31) of Up, making use of the expression (38) of y/p, and taking 
advantage of the identities ^1+^2 = i(io ~ ^3) ^or po,3 = and ti — t2 = i(^o + ^3) for pi 2 — (which can be 
established with the help of (39)), one finds that (jp is diagonal in the standard basis and has diagonal elements 

{n\ap\^^) = %t^(1 + (-l)V') , (uVplu) = ^^y^i^ + (-1)^) (53) 

with i,j = 0,l,i^ j. Here, r =- ±2Re{(01|$)($|10)} if P0P3 = and r = otherwise, and r' = ±2Re{(ll|<E>')(^'|00)} 
if P1P2 — and r' == otherwise. If poP3 = (respectively p\P2 = 0) then r (respectively r') can take arbitrary values 
between —1 and 1, since 1$) (|$')) is an arbitrary normalized vector in the two-dimensional space given in (52). 

The case |ci| > |c2|,|c3| (respectively |c2| > |ci|,|c3|) can be deduced from the previous case by the following 
symmetry argument. One can transform p into a state p' of the form (36) with a vector c' = (c3,C2,Ci) (respectively 
c' ~ (ci,C3,C2)) satisfying jcgj > |c']^|,|c2| by means of the unitary conjugation p' = UpW with U — e^^'^^ (E) e'^'^^ 
(respectively U — e~'~'^i e~'~'^^). By invariance of the Bures distance under unitary conjugations (property (ii) 
of Sec. II), such a transformation does not change the maximal fidelity, so that Fa{p) — Fa{p') = (1 + 63)72 = 
(1 + 6m^^^)/2 with TOmax = 1 (respectively rrimax = 2), in agreement with (45). Moreover, according to (51) the states 
CTp are related to the closest states a pi to p' by ap = U*<TpiU. Since api has diagonal matrix elements in the standard 
basis, given by (53) upon the replacement of ^3 by ig ~ immaxi this yields 

(^0 + immax) 



pr X ^(l + r')|ao,/3o)(ao,/3o| + (l-r')|ai,/3i)(ai,/3i 

+ ^*''~p17T^ {{l + r)\ao,l3,){ao,l3^\ + (1 - r)|ai,/3o)(ai,/3o|) (54) 

where we have set |ai,/3j) — \ai) (E) \/3j) with \ai) — \/3i) the eigenvectors of the Pauli matrix cr„i^^^, and 

{{(0, 0)} if all p^ are nonzero 

[-1,1] X{0} if poPm„ax = and p„ > 0, TO 7^ TOmax (55) 

{0} X [-1, 1] if poPm„ax > and P1P2P3 = 0. 

The explicit expressions of the eigenvectors of <Jm^^^ in the standard basis read 

\i) if Icsl > |ci|,|c2| 

|a.) = |A) = <( 2-1/2(10) + (-l)''|l)) if|ci|>|c2|,|c3| (56) 

^2-1/2(10) +i(-l)^|l)) if|c2|>|ci|,|c3|. 

One checks on formulas (54) and (56) that when c does not pertain to one of the faces of the tetrahedron (so that 
P0P1P2P3 > 0) and is such that \cm\ is maximum for a single component to, — rrimax, then ap is unique, has maximally 
disordered marginals, and is of the form 



4V 2Fa{p) 



^p = 7 ( f ® f ^ 77rT3 '^'"max «) cTm^ax j (57) 



(we have used here the identities Fa{p) = 4(tQ + t^„ ) and 16ioim = o,m + Cm, which follow from (37), (39), and 
(41)). The maximal fidelity Fa(p) in the denominator in (57) is given by (45). Thanks to (42), one gets 



ai 



I (V/(1+Cl)2-(C2-C3)2 - V/(1-Cl)2-(C2+C3)2) (58) 



16 

with analogous formulas for 02 and 03 obtained via a permutation of indices. 

We now turn to states p with vectors c such that |cm| is maximum for exactly two components c.^. For instance, 
let us assume ci = ±C2 and |ci| = |c2| > |c3| (the other cases will easily follow by a state rotation as above). Then 
61 = 62 > 63, see (47). Hence any vector u°p'((/)) = cos(t>ex + sa\(l>ey in the (a:;Oy)-plane defines an optimal direction 
(see Sec. VB). In other words, infinitely many basis {|q;°^ )} maximize the optimal success probability. The situation 
is similar to the case of pure states with degenerate Schmidt coefficients discussed in Theorem 1. The non- uniqueness 
of u°P* comes from the symmetry of p. For indeed, one finds from (36) that p = U±{(f)pU±{~(f) for any (j), with 
U±{4>) = e~''2'^3 (g) eT'2-<^3. By using the identity e'^a'^^^-^g-i-j-'^a — ^^, -^here u' is related to u by a rotation around 
the z-axis with the angle —(f>, it follows that 

A(u°P*(0)) = C/±((/<)A(e,)C/± (-(/)) . (59) 

As a result, the spectral projectors 11°^ (cj)) of A (u°p' (</))) and thus the closest A-classical state (31) corresponding to 
u°P*((^) are obtained from the corresponding projectors and A-classical states for u°p* = e^; by a unitary conjugation 
by U±{(l>). We conclude that all the A-classical states ap satisfying (51) are given 

ap{(t)) = -(10 1+ {cos(j)cri +sm<j>a2) <E) (cos(/)(Ti ±sin(/)(T2)j (60) 

if c does not belong to a face of the tetrahedron. More generally, the closest states a-p are of the form (54) with 

r2-i(|0) + (-l)V<^|l)) r2-^(|0) + (-l)^e±'-^|l)) ifci=±C2, |ci|>|c3| 

|a,) = <^e-'i-^|z) ,m^\e^'^^^i) if ci = icg, |ci| > |c2| (61) 

[e'^^^i) [e±'i-MO ifc2=±C3, |c2|>|ci| 

for arbitrary angles < (f) < 2tt and < 9 < 27r. 

Finally, let us study the states p with |ci| = |c2| = Icsj. This is the case for instance for the Werner states. Let 
^2,^3 G {~ljl} be defined by C2 = e2Ci and C3 — e^ci. Then p is invariant under rotations with arbitrary angles 9 

and 0, 

p=Ue,,eAd,^)pUe,.eAO,<l>)^ , t/e„e3 (0, 0) = e"'^-^ e"'^^^ ® 6"^^ ^'^^ g-i.3 f ^2 _ (62) 

Moreover, the eigenvalues of A(u) are independent of u, since &i = &2 = ^3, see (43). Thus the optimal direction u°p* 
is completely arbitrary and any orthonormal basis {|a°^ )} of C^ is optimal. The closest A-classical states to p are of 
the form (54) with 

\a,) = e-'^'^^e-'^'^^lz) , |ft) = e-'^^^'"^e-'^^^'"^|i) . (63) 

The expressions in the standard basis of \ai) and \/3i) in (61) and (63) can be obtained from the relations 

g-i.2|<T3g-ie3f <T2 |o) ^ COS - |0) + 63 e'""*' sin - |1) 

g-ie2|<T3g-ie3f<T2|^^ ^ sin- |0) - €36''^^'^ COS- 1 1) (64) 

where we have disregarded irrelevant phase factors. 

It should be noted that the closest ^-classical states (54) are in fact classical states. This means that the distance 
of a state p with maximally disordered marginals to Ca coincides with its distance to the set C of classical states, i.e.. 



D{p) = Da{p) = Db{p) = dsip^ap) = 2(1 - ,Jl±^^^ 



(65) 



as could be expected from the symmetry of p under the exchange of the two qubits. 

In Fig. 1 (b),(c), some examples of states p in the tetrahedron and their closest states ap are represented. Outside 
the dashed lines on which two or more components c^ have equal moduli, ap is specified by a vector s lying on the 
closest coordinate semi-axis to c for the usual distance in R"^ (but s is not the closest vector to c on that semi-axis 
for this distance). 

The closest yl-classical states for the Hilbert-Schmidt distance ^2 to a state p of the form (36) can be found by 
using the results of Ref. [21]. They are given by 



3 

4 

l^m—l 



1 ^ 

'f'> = -(1 (E) 1 + ^ UiCmUmm «) (Jmj (66) 
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with an arbitrary unit vector u G R^. In contrast with the Bures distance, all states p G T have infinitely many closest 

(2) (2 

Tp . Moreover, there are three states Up 



(2) (2) 

A-classical states ap . Moreover, there are three states Up of the form (36) specified by vectors s^ = Cmem lying on 
the coordinate axes tti = 1, 2, 3, irrespective of the order of the moduli \cm\- Since tr[((Tp )^] = (1 + X]m '^m''^m)/4 ^ 



(2) 

1/2, the closest states ap are always mixed. In particular, unlike for the distance ds (see Theorem 1), the closest 
A-classical states to the maximally entangled pure states |$^) or |^^) are not pure product states. 

VI. CONCLUSIONS 

We have established in this paper a link between the ambiguous quantum state discrimination task and the problem 
of finding the minimal Bures distance of a state p to a state with zero quantum discord. More precisely, the maximal 
fidelity between p and an A-classical (i.e., zero-discord) state coincides with the maximal success probability in 
discriminating ua states p°^ with prior probabilities 7]°^ , given by Eq.(18) and depending on an optimal orthonormal 
basis {|a°^ )}"=i of subsystem A (Theorem 3). The closest A-classical states to p are, in turn, given in terms of this 
optimal basis and of the optimal von Neumann measurements in the discrimination of {p°^ , 7]°^ } (Theorem 4) . By 
solving the quantum state discrimination problem we have found explicit formulas for the Bures-GQD and for the 
closest A-classical states to a given two-qubit state with maximally disordered marginals. 

Since the quantum state discrimination task can be solved for ua = 2 states, one can in principle by our method 
determine the GQD and the closest A-classical states for arbitrary two-qubit states, or even for arbitrary states of a 
bipartite system composed of a qubit A and a subsystem B with Hilbert space dimension ub > 2. An open question 
for two qubits concerns the bound (50) between the quantum discord and the Bures-GQD, which holds for two-qubit 
states with maximally disordered marginals but could be more general. For subsystems A with higher Hilbert space 
dimensions ua > 2, several open issues deserve further studies. Firstly, it would be desirable to characterize the 
"most quantum" states having the smallest fidelity to A-classical states when ua > ub- Secondly, it is not excluded 
that the specific quantum state discrimination task associated to the minimal Bures distance admits an explicit 
solution. Thirdly, the relation of the GQD with the geometric measure of entanglement should be investigated; in 
particular, there may exist some inequality analog to the monogamy relation [12] between the quantum discord and 
the entanglement of formation. 

Our results on the closest A-classical states may shed new light on dissipative dynamical processes involving deco- 
herence, i.e., evolutions towards classical states. In fact, our analysis allows to determine the geodesies passing through 
a given state po with a non-zero discord. Such a geodesic contains all states p having the same closest A-classical 
state ffpg as pQ. It would be of interest to compare these Bures geodesies with the actual paths followed by the density 
matrix in specific physical examples. 
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Appendix A: A necessary and sufficient condition for the optimal success probability to be equal to the 

inverse number of states 

Let {pi}"^i be a family of ua states on "H with prior probabilities rji, where ua = N/ns is a divisor of dim('H) = N. 
We assume that the pi have ranks rank(/9i) < ns for any i. Let Pg^ ^' '({pi,??i}) be the optimal success probability 
defined by Eq.(17). We prove in this appendix that Pg^ ^' '{{pi, r/i}) — I/ua if and only if r/i — I/ua for any i and 
all states pi are identical. 

The conditions rji — I/ua and Pi — p are clearly sufficient to have Pg^ ^' '{{pi, rji\) — I/ua (a measurement cannot 
distinguish the identical states pi and thus cannot do better than a random choice with equal probabilities). We need 
to show that they are also necessary conditions. Let us assume Pg^ ^' ■({pi,??i}) ^ l/"-A- The equality of the prior 
probabilities rji is obvious from the bounds Pg^ ^' '({pi,'7i}) > Vmax and ?7niax > ^/tt-a (see Sec. IV A, which implies 
Vmax — l/^A ~ Vi for ^-H «• As a result and by using our hypothesis, any orthogonal family {ni}"J'j^ of projectors 
of rank ub satisfies '^^^^tr (Hi pi) — Ti^Pgdp^, 77^}) = 1. We now argue that the states pi have ranges contained in a 
common ns-dimensional subspace. In fact, let Vi be the ns-dimensional subspace of "H containing the eigenvectors of 
pi associated to the ub highest eigenvalues, and let us denote by Hi the projector onto Vi. Then ranpi C Vi (since 
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we have assumed rank(pi) < ub) and thus pi — Hipi- This imphes 1 = '^^tr{IliPi) > tr{Ilipi) = 1. It follows that 
tr(n2P2) — for any projector 112 of rank ns orthogonal to Hi. Hence pi and p2, and similarly all p^, i = 3, . . . , n^^, 
have ranges contained in Vi. In order to show that all the states pi are equal, we further consider, for each 1 < A: < ub, 

(k) 

some ns-dimensional subspace Vl containing the eigenvectors associated to the k highest eigenvalues Ai, . . . , A^ of 
Pi, the other eigenvectors being orthogonal to V} (then v}"^ = Vi). We also choose a n^-dimensional subspace 
W[''^ C n orthogonal to V^") such that wf ^ ® v[''^ D Vi. Let {nf ^}"4 be an orthogonal family of projectors of 
rank ub such that Hi ' and 112 ^^^ the projectors onto V} and W} , respectively. Then 

1 = Y, tr(nf V») = tr(nJ'=Vi) + tr[(l - U^t^)p2] = 1 + Ai + . . . + Afe - tr(n('V2) (Al) 

i 

where we used ^^ 11^ = 1 and ran(pi) C W} ® V} in the second equality. By virtue of the min-max theorem, 

tr(n]^ P2) is smaller than the sum of the k highest eigenvalues of p2- The identities (Al) imply that this sum is larger 
than the sum Ai + . . . + A^ of the k highest eigenvalues of pi . By exchanging the roles of pi and p2 we obtain the reverse 
equality. Since moreover k is arbitrary between 1 and ub, it follows that pi and p2 have identical eigenvalues. By 

using (Al) again, tr(n-j^ P2) is equal to the sum of the k highest eigenvalues of p2, implying that the k corresponding 

(k) 
eigenvectors of p2 are contained in Vj^ . Since k is arbitrary, this proves that pi and p2 have identical eigenspaces. 

Therefore pi = p2. Repeating the same argument for the other states p^, i > 3, we obtain pi = . . . = Pn^- 
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